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1. Introduction 



The inhomogeneous Bethe-Salpeter (BS) equation in Minkowski space §\§ provides a covari- 
ant four-dimensional description of two-body scattering states. In the case of scalar particles it has 
the form: 

d 4 p' 



F(p,p s ;P)=K(p,p s ;P)-iJ ^ 



K(p,p';P)F(p', Ps ;P) 



P + p') -m 2 + ie (\P-p') -m 2 + ie 



(1.1) 



The kernel K for the one-boson exchange model is given by: 

„2 



k(p,p';p) 



(p - p') 2 - H 2 + ie 
We introduce the coupling constant a related to g by: 

g 2 = \6nm 2 a 

and use the partial wave decomposition according to |§] : 

oo 

F(0) = 16n£(2Z + l)F/fl(cose) 

1=0 



(1.2) 



In the center of mass frame, P = 0, Po = yfs = 2e Ps = lym^+p 2 and for a given incident momen- 
tum p s , the partial wave off-mass shell amplitude F\ depends on two scalar variables po and \p\. It 
will be hereafter denoted by Fi(po,p; p s ) setting p = \p\, p s = \p s \. We consider for simplicity the 
S-wave equation which reads: 

4i 



4i r°° 2 r°° 

Fo(po,p;ps)=Ko(po,p;ps) — 2 P d P d P0 

Ko(po,p;p ,p')F {p'o,p''>Ps) 



(1.3) 



(p'l + 2p' Q e Ps + p 2 - p' 2 + ie) (p'l - 2 Po e Ps + p 2 - p' 2 + ie 



with 



k (po,p;p'o,p' 



1 f 1 



32tt (p - p' Q ) 2 - (p 2 - Ipp'z + p' 2 )-p. 2 + ie 

2 |„ , 11 ;«„„2 



am', + ictizm 



■lOE 



+ 



4p// | T7 — 1 1 4pp 

C/(7]) = 



(1.4) 



1, if |tj| < 1 
0, if 1 77 1 > 1 



and 



(Po - Po) 2 - P 2 - P' 2 - H 
Ipp' 



The on-shell amplitude F,°" = Fi(po = 0,p = p s ; p s ) determines the phase shift according to: 

(1.5) 



2i 
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The knowledge of this function in the entire domain of its arguments - i.e. the off-shell amplitude - 
is mandatory for some interesting physical applications, like for instance computing the transition 
e.m. form factor y*d — > np or solving the BS-Faddeev equations. As far as we aware, this quantity 
has not been obtained until now. 

The numerical solution of the BS equation in Minkowski space is complicated by the existence 
of singularities in the amplitude as well as in the integrand of These singularities are inte- 

grable in the mathematical sense, due to is in the denominators of propagators, but their integration 
is a quite delicate task and requires the use of appropriate analytical as well as numerical methods. 

To avoid these singularities, the BS equation was first solved in Euclidean space. These so- 
lutions provided on-shell quantities like binding energies and phase shifts [g]. However we have 
shown [Q] that the Euclidean BS amplitude cannot be used to calculate electromagnetic form fac- 
tors, since the corresponding integral does not allow the Wick rotation. One therefore needs the BS 
amplitude in Minkowski space. 

This amplitude has been computed for a separable kernel (see []|] and references therein). For 
a kernel given by field theory rules - ladder and cross ladder - the Minkowski BS amplitude was 
first obtained in our preceding works []6| ^ in the bound state problem. To this aim, we developed 
an original method based on the Nakanishi integral representation of the amplitude. A similar 
method for the scattering states has been proposed in []8|] although the numerical solutions are not 
yet available. 

We present in this paper some details of a new method [^] providing a direct solution of the 
original BS equation. It is based on a carefull treatment of the singularities and allows to compute 
the corresponding off-shell scattering amplitude in Minkowski space. We will give the low energy 
parameters in the case of spinless particles and ladder kernel. 

2. Method 

There are four sources of singularities in the r.h.-side of the BS equation Ql.l[ ) which are 
detailed below. 



(i) The constituent propagators in (1.3) vs. p' Q have two poles, each of them is represented as: 



p' — a — ie P'q — ci 



PV^— + iKd(p' -a) 



where PV means the principal value. In the product of four pole terms, the only non vanishing 
contributions come from the product of four PVs without delta-functions, from the terms with 
three PVs and one delta-function and from the term with two PVs and two delta's. After partial 
wave decomposition the 4D integral BS equation ( |Q| ) is reduced to a 2D one, eq. (L3). Integrating 



in ( |l.3| ) over p' , we obtain in addition to the 2D part, a ID integral over p' and a non integrated 
term. The singularities due to the PVs are eliminated by subtractions according to the identity: 



PV 



fiPoWo 
Pq -a 1 




The integrand in r.h.-side is not singular. 
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( ii) The propagator of the exchanged particle has the pole singularities which, after partial wave 
decomposition, turn into logarithmic ones, eq. ( |l.4[ ). Their positions are found analytically and the 
numerical integration over p' variable is split into intervals between two consecutive singularities, 
namely: 



/ [...}d P ' Q = / [...}d P ' + / [...}dp' Q + .. 

JO Jo Jsing, 



'sing i 

Each of these integrals is made regular with an appropriate change of variable. We proceed in a 
similar way for the p' integration. 

(Hi) The inhomogeneous (Born) term is given by the ladder kernel and is also singular in both 
variables. The positions of these singularities are analytically known. 

(iv) The amplitude Fo itself has many singularities, among which the strongest ones are re- 
sulted from the Born term Ko(po,p; p s ). This makes difficult its representation on a basis of regular 
functions as well as its numerical integration. To circumvent this difficulties we made the re- 
placement Fo(po,p; p s ) = Ko(po,p; p s )fo{po,P\ Ps), where fo is a smooth function. After that, the 
singularities of the inhomogeneous term are canceled. We obtain in this way a non-singular equa- 
tion for fo which we solve by standard methods. Then we restore the BS off-mass shell amplitude 
Fq in Minkowski space. 



3. Numerical results 

We first applied this method to solve the bound state problem by dropping the inhomogeneous 
term in (|L~l|). The binding energies coincide, within four-digit accuracy, with the ones calculated 
in our previous work |Q] and with the Euclidean space results. 

Solving eq. ( |l.3| ), the S-wave off-shell scattering amplitude Fo is calculated and the phase shifts 
are extracted by means of eq. (L5). Above the first inelastic threshold , pJ(jU) = \Jm\i + pl 2 /4, 



this phase shifts have an imaginary part, which is also found. By performing a Wick rotation in 
( |L~3| ) and taking into account the contributions of singularities - which, in contrast to the bound 
state case, are crossed by the rotation contour - we derived an Euclidean space equation similar to 
one obtained in [^]. The phase shifts found by these two methods - i.e., solving eq. ( 1.3) and the 



Euclidean space equation - coincide with each other within 3-4 digits. Furthermore, the imaginary 
part of the phase shifts vanishes with high accuracy below threshold. The unitarity condition is 
not automatically fulfilled in our approach, but appears as a consequence of handling the correct 
solution. It thus provides a stringent test of the numerical method. Our results reproduce the phase 
shifts given in [^] within the accuracy allowed by extracting numerical values from published 
figures. 

Figure[l] P anei ) shows the phase shifts calculated via BS equation (solid curve) and via 
the Schrodinger one with the Yukawa potential (dashed curve) for the constituent mass m = 1, 
exchange mass pL = 0.5 and coupling constant a = 0.5. Right panel shows the same results for 
a = 1.2. For this value of a, there exists a bound state. Therefore, according to the Levinson 
theorem, the phase shift starts at 180 degrees. One can see that the difference between relativistic 
and non-relativistic results is considerably large, specially for small incident momentum. This 
difference increases with a. 
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Figure 1: Left panel: phase shift calculated via BS equation (solid curve) compared to non-relativistic 
results (dashed curve) for jj, = 0.5 and a = 0.5. Right panel: the same as for the left panel but for a = 1 .2 




Figure 2: Left panel: imaginary part of the phase shift (solid) and squared modulus of the S-matrix (dashed). 
Right panel: scattering length oq vs. a for fi = 0.5 calculated via BS (solid) and Schrodinger equations 
(dashed). 



Fig. U (left panel) shows the imaginary part of the phase shift (for a = 1.2) which auto- 
matically appears when the incident momentum exceeds the threshold value for creation of the 
exchange meson. For m = \, \l = 0.5 this value p thresh = 0.75. Simultaneously, the modulas of 
two-body S-matrix differs from 1. For p s = 1.118 the second threshold, for creation of two mesons, 
is open. It also contributes in this curve. 

Fig. U (right panel) shows the scattering length c?o as a function of the coupling constant 
a obtained with BS and non-relativsitc Schrodinger equations. In the vicinity of a ps 0.8 (for 
Schrodinger) and awl (for BS) the coupling constant crosses the critical value corresponding 
to the appearance of a bound state. At this point the scattering length becomes infinite and then 
changes the sign. 

We have displayed in Fig. || the real (left panel) and imaginary (right panel) parts of the off- 
shell scattering amplitude Fo(po,p; p s ) vs. po and p calculated for p s = /J. = 0.5. Its real part shows 
a non trivial structure with a ridge and a gap resulting from the singularities of the inhomogeneous 
term. Its on-shell value Fq" = Fq(pq = 0,p = p s ;p s ), determining the phase shift calculated pre- 



5 



Scattering states in Bethe-Salpeter equation 



V.A. Karmanov 



viously, corresponds to a single point po = 0,p = p s on these surfaces. Our calculation, shown in 
Fig. ^j, provides the full amplitude Fo(po,p;p s ) in a two-dimensional domain. It cannot be found 
from the Euclidean equation. Computing this quantity is the main result of this work. 



Figure 3: Real (left) and imaginary (right) parts of the off-shell amplitude F(po,p; p s ) for p s = 0.5, p, = 0.5. 

4. Conclusion 

We solved the BS equation for the scattering states in Minkowski space for the ladder kernel. 
The off-mass-shell amplitude is found for the first time. Coming on mass shell, we obtain the 
phase shifts which coincide with ones calculated by other methods. They considerably differ, even 
at low energy, from the non-relativisttic phase shifts calculated by the Schrodinger equation. Above 
the meson creation threshold the inelasticity appears which is also calculated. The off-mass-shell 
amplitude can be used to calculate the transition form factor and as an input in the three-body 
BS-Faddeev equations. 
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